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PULSE PHASE RESOLVED ANALYSIS OF THE HIGH MASS X-RAY BINARY CENTAURUS X-3 OVER TWO BINARY 

ORBITS 

Slawomir Suchy', Katja Pottschmidt''- \ Jorn Wilms'*, Ingo Kreykenbohm''^, Gabriele Schonherr', 
Peter Kretschmar**, Vanessa McBride', Isabel Caballero^, Richard E. Rothschild', Victoria Grinberg'" 

ABSTRACT 

We present a detailed analysis of observations of the high mass X-ray binary Cen X-3 spanning two consec- 
utive binary orbits performed with the RXTE satellite in early March 1997. During this time Cen X-3 had a 
luminosity of L2-iokeV ~ 4-5 x lO^'^ergs"' and a pulse period of 4.814 s. The PCA and HEXTE light curves 
both show a clear reduction in count rate after mid-orbit for both binary revolutions. We therefore analyze two 
broad band spectra for each orbit, before and after mid-orbit. Consistent with earlier observations these four 
joint PCA and HEXTE spectra can be well described using a phenomenological pulsar continuum model, in- 
cluding an iron emission line and a cyclotron resonance scattering feature. While no strong spectral variations 
were detected, the second half of orbit 2 shows a tendency toward being softer and more strongly absorbed. 
In order to follow the orbital phase-dependent evolution of the spectrum in greater detail, we model spectra 
for shorter exposures, confirming that most spectral parameters show either a gradual or sudden change for the 
second half of the second orbit. A comparison with a simple wind model indicates the existence of an accretion 
wake in this system. We also present and discuss high resolution pulse profiles for several different energy 
bands, as well as their hardness ratios. PCA and HEXTE spectra were created for 24 phase bins and fitted 
using the same model as in the phase averaged case. Systematic pulse phase-dependent variations of several 
continuum and cyclotron line parameters were detected, most notably a significant increase of the cyclotron 
line energy during the early rise of the main peak, followed by a gradual decrease. We show that applying a 
simple dipole model for the magnetic field is not sufficient to describe our data. 

Subject headings: X-rays: stars — X-rays: binaries — stars: pulsars: individual (Cen X-3) — stars: magnetic 
fields 



1. INTRODUCTION 

The pulsar nature of the accreting high mass X-ray bi- 
nary (HMXB) Centaurus X-3 was first dis covered by UHURU 
dSchreier et al.ll 19721 iGiacconi et al.|[T971D . measuring a spin 
period of ~ 4.8 s and an orbital period of 2.1 days. The bi- 
nary system consists of a neutron star with a mass of 1.21 ± 
0.21 M0 accompanied by an O 6-8 III supergian t star with 
a ma ss of 20.5 ± 0.7Mq jHutchings et al'l 119791: lAsh et all 
119991) . The distance of the binary system has been estimated 
to be roughly 8 kpc with a lower limit o f 6.2 kpc (jK rzeminskj 
1974 and an ec centricitv of < 0.0016 C Bildsten et al.lll997l) . 
Day & Tennanj (Il991i) used dust scattering measurements to 
derive a distance of 5.4 ±0.3 kpc, a result only marginally 
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consistent with' Krzeminskil (1197 4') within given uncertainties. 
In order to allow for easier comparison with earlier publica- 
tions we will adopt the widely used distance of 8 kpc. HMXBs 
generally show a strong stellar wind which can supply a con- 
tinuous tidal stream to the compact object when the com- 
panion star is close to filling its Roche lobe (Petterson 1978t 
Stevens 1988; Blondin, Stevens & Kallman 1991). Dav et aQ 
7r993l) discussed that in a system like Cen X-3 an X-ray 
excited wind would be sufficient as a possible mechanism 
for mass transfer able to sustain a constant mass flux. A 
system with a combined mas s transfer mechanism is, e.g., 
GX 301-2 (Leahy 2002; Leahy.&Kostka 2007), where the 
folded RXTE /ASM light curve was best modeled by a stellar 
wind with an additional gas stream as described by ,Stevensl 
(Il988l) . In additi on evidence for an accretion disk in Cen X- 
3 can be found: iTiemkes. van Paradiis & ZuiderwiikI ( Il986h 
compared the observed light curve of Cen X-3 with mod- 
eled binary light curves and achieved best results when in- 
cluding an accretion disk. The observed overall spin-up trend 
of 1.135 ms/yr with fluctuations on ti mescales of a few years 
(iTsunemi. Kitamoto & Tamurall 19961) also indicates the pres- 
ence of an accretion disk, attributed to Roche lobe overflow. 
A realistic picture of the mass transfer mechanism in Cen X-3 
would therefore be a combination of disk ac c retion with an 
excited stellar wind, as proposed by lPettersonI (Il978h . 

An absorption feature in the X-ray spectrum at ~30keV, in- 
terpreted as a cyclotron resonance scattering feature (CRSF) 
also known as cyclo tron line, w as first observed in Cen X- 
3 with Ginga by Nagase et alj ( [1992) . BeppoSAX and 
i?XT£/HEX TE confirmed this featur e at 28keV and 30keV 
respe ctively dSantangelo et al.l Il998t iHeindl & Chakrabartvl 
119991) . It is assumed to be the fundamental line and so far 
no higher line harmonics have been observed for this source. 



2 



Suchy et al. 



CRSFs are generated due to the interaction of photons with 
electrons which are quantized in Landau levels, the discrete 
energy states populated by electrons in strong magnetic fields. 
Photons of energies corresponding to the separation of the 
Landau levels are absorbed and re-emitted by the electrons. 
The net effect is that photons from the accretion column above 
the magnetic poles that have energies near the CRSF cen- 
troid energy are scattered out of the line of sight and/or re- 
distributed in energy. These features can be used to measure 
directly the magnetic field strength in the line-forming region 
since the fundamental cyclotron line energy is given by 
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where B\2 is the magnetic field strength near the neutron star 
surfa ce, in units of 10'^ Gauss and z is the gravitational red- 
shift (ICanuto & Venturalll977h . Assuming a typical neutron 
star mass of 1.4M0 and neutron star radius of 10 km gives 
z = 0.3. Cen X-3's cyclotron line at ^ 30keV then corre- 
sponds to a mag netic field of ~ 3.4 x 10'^ G. 

Kohmura. Kitamoto & Toriil (12001) analyzed the PCA data 
of this observation in the energy band 2- 13keV which 
has a high time resolution of 15.625 ms. They saw a de- 
lay in the time variability of the iron K band by 0.39 ± 
0.10 ms, and concluded that these photons were reprocessed 
away from the neutron star and are probably part of the 
matter accreting onto the N S. A preliminary analysis by 
iHeindl & Chakrabartvl (IT999I) of the HEXTE data presented 
in this work concentrated on basic pulse phase resolved spec- 
troscopy, dividing the pulse profile into four broad phase 
bins, and provided a first indication of a possible increase 
of the cyclotron line energy during the rise of the main 
pulse. Systematic variations of the CRSF centroid en- 
ergy over pulse phas e have also been observed in othe r 
sources, e.g., HerX-1 (Soong et al. 1990; Gruber et al."2001'). 
Vela X-1 (La Barbera et al. 2003; Krevkenbohm e t al. 2002J 
4U 03 52+309 ( Coburnet al. 2001), 4U 153 8-52jClark et al.1 
[T990h . and GX 301-2 dKreykenbohm etani2004l) . The ori- 
gin of these changes is not yet well understood. Straightfor- 
ward mechanisms like the angular dependence of relativistic 
corrections to the cyclotron energy may not be sufficient to 
explain the observed variations. Higher harmonics rarely, if 
ever, are integer multiples of the fundamental line energy. In 
some sources this can be explained by invoking the proper 
non-integer line spacing expected from the relativistic cross 
section (e.g., Pottschmidt et al. 2005). However, more com- 
plicated mechanisms such as a B-fiel d gradient at the mag- 
netic pole, have also been invoked (iNishimural l2005h . In 
addition a correlation between CRSF energy and luminos- 
ity can be observed in several sources (e.g., Mihara '1995'; 
iMowlavi et al. 2006; L a Barbera e t al. 2003; Staubertetal 
I2OO7I: iNakaiima et al.ll2006h . 

In the following we further characterize the phase depen- 
dence of the line parameters in Cen X-3 using the excep- 
tionally long - spanning '^4days, i.e., two bina ry revolutions 
- RX TE observation of the source in 1997. iBurderi et alj 
(|2000), in the following referred to as BOO, observed Cen X- 
3 simultaneously with BeppoSAX overlapping 1 1 hours of 
RXTE observations during the first binary orbit. We will 
compare our results with BOO whenever possible. 

We use combined PCA and HEXTE data for the analysis 
which allows us to constrain the continuum shape of the spec- 
trum and the CRSF parameters much more strongly than with 
HEXTE alone. We compare our analysis with other cyclotron 



line sources and apply a newly developed physical model for 
the CRSF to the data. For the analysis over the orbit we com- 
pare with a simplified wind model to discuss the evolution of 
A^H through the observation. 

The remainder of this paper is organized as follows: ba- 
sic observation parameters as well as the data reduction are 
described in ^2] In ^ we describe the PCA and HEXTE 
light curves observed over the two binary orbits. Phase aver- 
aged spectra for the first and second halves of each orbit and 
their best fit models are presented in ^ as well as the tem- 
poral evolution of the spectral parameters with a resolution of 
roughly one RXTE orbit throughout the two Cen X-3 binary 
orbits, ^is dedicated to the pulse phase resolved analysis, 
including the presentation of pulse profiles in multiple energy 
bands and modeling of phase resolved spectra. The results are 
discussed in |6] 

2. OBSERVATION AND DATA REDUCTION 

RXTE obs erved Cen X-3 wit h the Proportional Counter 
Array (PCA; Jahoda et al." 2006) and the High Energy X- 
ray Timing Experiment (HEXTE; Rothschild et al. 1998) be- 
tween 1997 February 28 and March 3. The observation, pro- 
posal ID P20104, covered two consecutive binary orbits be- 
tween eclipses, excluding the time of eclipse between orbits. 
The total exposure time of 3 1 ks was split into 1 2 shorter data 
sets ("Obslds") of ^ 8 hours each. 

The PCA consists of five Proportional Counter Units 
(PCUs) covering a nominal energy range from 2 to 60keV. 
For ^ 90% of the time all PCUs were collecting data during 
the observation, with PCU 3 being turned off during the re- 
maining time. The PCA background was modeled using the 
"Sky-VLE" model. The HEXTE instrument consists of two 
independent clusters, each containing four NaI(Tl)/CsI(Na) 
phoswich scintillation counters, covering an energy range 
from 15 to 250 keV. Each cluster pointed alternately at the 
source and at the background, switching every 16 s. At any 
given time one of the clusters was on-source. 

The data were reduced using HEASOFT version 6.0. 6., ap- 
plying standard extraction criteria: a pointing offset of < 
0.01° from the nominal source position, an exclusion time for 
the South Atlantic Anomaly (SAA) of 30min, and a source 
elevation of > 10°. We imposed a maximum "electron ra- 
tio" of 0. 15 in order to filter for background flares, taking into 
account a slightly elevated electron ratio level due to the in- 
fluence of the bright target source. Depending on the anal- 
ysis method, further grouping and filtering of the resulting 
good time intervals (GTIs) were applied in order to derive the 
PCA and HEXTE light curves and spectra. Details are de- 
scribed below, including the respective choices of data modes 
and binning options. All HEXTE data products were dead- 
time corrected using the HEASOFT ftool hxtdead. Since 
precise absolute flux measurements are beyond the scope of 
this paper, correction for the PCA deadtime of a few percent 
has not been carried out. Systematic uncertainties of 0.5% 
were included in the PCA spectral data. All high time reso- 
lution products were converted into the frame of reference at 
baryce nter The orbital ephemeris presented by lBildsten et alJ 
( Il997h was adopted for the binary star corre ction. For spec tral 
modeling we used XSPEC version 11.3.2p (lArnaudl[T996l) . 

3. BROAD BAND ORBITAL LIGHT CURVE 

We extracted PCA and HEXTE light curves for all available 
ObsIDs. For the PCA, "standard2f" mode data were used, se- 
lecting the well calibrated energy range from 3-17 keV and 
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Fig. 1 . — (a) Background subtracted 3-17 keV PCA (top) and 17-100 keV HEXTE (bottom) light curves (bintime 128 s). Shaded region indicates the eclipse 
of the compact object. Error bars are not plotte d bu t are smaller than any variations visible. Horizontal lines above the data indicate the time intervals used for 
creating the time averaged spectra analyzed in j|4. II (b)-(e) Evolution of selected spectral parameters over two binary orbits. See text for model details. While 
the spectrum is stable - away from enhanced (pre-)eclipse - over the first orbit, the behavior changes significantly after the mid-orbit break of the second 
orbit. Data points from spectral analysis during eclipse are ignored. 
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the highest possible time resolution for this data mode of 
16 s. Corresponding background model light curves were cre- 
ated. In the case of HEXTE we use science event mode data 
from 17-100 keV, with the source not being significantly de- 
tected at higher energies. Since the HEXTE source and back- 
ground measurements for each cluster are performed quasi- 
simultaneously due to the 16 s rocking cycle, the high time 
resolution background for HEXTE light curves with less than 
16 s time resolution is based on interpolating the measured 
background rates. For each 1 s time-resolved data bin in a 
16 s segment we performed a linear interpolation of the back- 
ground, using the neighboring 16s background pointings". 
As HEXTEs cluster B only has 3 fully functional detectors, 
the count rate for this cluster was normalized to cluster A. 

Fig.lTJj shows the background subtracted PCA and HEXTE 
light curves (both rebinned to a resolution of 128 s for display 
purposes). The PCA light curve has been normalized to one 
PCU and the HEXTE light curve consists of the summed con- 
tributions of both clusters and has been multiplied by a factor 
of nine, for display purposes. The same general features are 
visible above and below 17 keV, indicating the absence of ex- 
treme spectral variability over the binary orbit. For orbital 
phase (/)orb 0.5 both orbits show a clear drop in the count 
rates. In addition the source was fainter overall during the sec- 
ond orbit, with the brighter first half of orbit 2 approximately 
at the same level as the less bright second half of orbit 1 . The 
eclipse covers ~ 20% of the orbit al period. Pre- e clipse dips 
can be obser v ed in both orbits. .Takeshima et al.l (|1991|) and 
iNagase et al.l (Il992h also observed a drop in luminosity for 
the second part of the orbit, calling it a pre-eclipse dip. In our 
case we distinguish between before and after mid-orbit break 
and use the term pre-eclipse dip for relative short count rate 
drops, slightly longer than one satellite orbit. 

In the first orbit one single pre-eclipse dip is visible around 
001-b = 0.75 (JD 2450508.4) where the PCA rate drops by a 
factor of two. Unfortunately, no HEXTE data are available 
in the archive for the egress of this pre-eclipse because the 
instrument was not turned on (PCA data are available). In 
the second orbit dipping activity starts around the same ^orb 
(JD 2450510.8), however, all pre-eclipse dips are partially ob- 
scured by data gaps due to SAA passages and earth occulta- 
tions, so that no detailed analysis is possible. 

Based on modeling the averaged spectrum for the first part 
of orbit 1 (lA in Fig. ffland Table [B we find a 2- lOkeV 
unabsorbed flux of 7.0:l:o7 ^ lO'^ergcm'^s"' with the PCA. 
Taking the deadtime influence into account increases the 
measured flux by a few percent. The BeppoSAX observa- 
tion analyzed by BOO was performed between 1997 Febru- 
ary 27 19:45 and February 28 11:00 (UT) and thus overlaps 
with parts of data set lA of the RXTE observation. BOO 
find a somewhat lower post eclipse egress 2-lOkeV flux of 
5.7 X 10~^ergcm~^s~'. The 18% discrepancy in flux is not 
surprising, due to the uncertainty in the cross calibration be- 
tween the PCA and the BeppoSAX. We point out that the 
fluxes derived from the HEXTE instrument are 14% lower 
than those derived with the PCA (see chexte in Table ^ 
and thereby agree much better with BeppoSAX. Based on 
the PCA flux measurement we derived a 2-lOkeV luminos- 
ity of 5.4;!;q g X 10^^ ergs"', applying a distance of 8 kpc (BOO: 
4.4x lO-'^ ergs"'). The fainter second orbit shows a 2-10 keV 



" More precisely the background is measured during 12 s, preceded and 
followe d by 2 s durin g which the cluster is moving between on- and off- 
source l IRothschild et al.ll998i) . 



TABLE 1 

Best fit parameters for phase averaged spectra over first 
and second halves of each binary orbit. 
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Note. — The data sets 1A-2B are defined in Fig.[T] Parameters listed 
are the hydrogen column density Nn; the Fermi-Dirac cutoff and folding 
energies i^cut 

and fifoij; the CRSF energy Ecyc ? width (Jt^yt; , and depth Tcyc 1 
the power law index T and strength Api; the Fe K emission line energy 
iJpe, width CTpe, and strength Ape; and the 13 keV emission feature energy 
£^i3kcV. width o-i3i;ev, and strength A 13 i^ev- Uncertainties are quoted at the 
90% confidence level assuming independent parameters (see text for exist- 
ing coiTelations). The flux normalization of the HEXTE instrument with 
respect to the PCA, chexte, and the luminosity in the 2-10 keV energy 
range derived from the PCA, Z,2-lOkeV> are also given. 
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during its first half, corre- 
lO-'^ergs"'. The lumi- 
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sponding to a luminosity of 4.3 
nosities after the mid-orbit drop are 4.4;|;['j x 10^^ ergs"' and 
4.0;!;[ g X 10^^ ergs"' for orbit 1 and 2, respectively, i.e. the rel- 
ative flux change in the second orbit is a factor of two smaller 
than in the first orbit. 

4. PHASE AVERAGED SPECTROSCOPY 

4.1. Time averaged spectral parameters 

For each available Obsid "standard2f" mode PCA spec- 
tra and standard science event mode HEXTE spectra were 
extracted, including the corresponding background spectra. 
All spectra within one of the four time intervals indicated in 
Fig.Hh were added. The resulting time averaged spectra are 
count rate selected in the sense that they are characterized by 
episodes of different average count rate levels within a given 
binary orbit, i.e., for each of the two orbits the first interval 
spans the bright phase up to the mid-orbit break, while the 
second interval spans the time up the on-set of pre-eclipse dip- 
ping. Time averaged broad band source spectra were modeled 
taking background subtracted PCA spectra from 3-23 keV 
and background subtracted HEXTE spectra from 17-100 keV 
into account, rebinning data above 60keV by a factor of five. 
As the CRSF centroid energy is at ^ 30keV and a contami- 
nation of our data due to the Xenon K-edge at 33keV would 
be possible, we decided not to increase the PCA range up to 
60keV as described by Rothschild et al. ( 2006), Appendix B. 

The basic continuum we use throughout the paper is an ab- 
sorbed power law with index F, modified by a Fermi-Dirac 
cutoff characterized by cutoff and folding energies E^ut and 
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Fig. 2. — (a) PCA and HEXTE phase averaged counts spectra, accumulated 
over the first part of the first binary orbit (lA in Fig.[T]and Table [T), and 
best fit model, (b) Residuals for an absorbed power law with a Fermi-Dirac 
cutoff, and a Gaussian Fe line, (c) Improved residuals obtained by adding 
a cyclotron line at ~ 30 keV. (d) Best fit residuals obtained by modeling the 
broad residual around 13 keV with an additional Gaussian emission line. The 
reduced values are 24, 3.4, and 0.97 from (b) to (d). See text for a detailed 
description of the model components and Table[T]for the best fit parameters. 

Efoid (local XSPEC model f dcut. lT^nakalll986h 

FiE) = AE-^- \ , (2) 

a phenomenological model co mmonly applied to describ e 
spectra of accreting pulsars, e.g..'Krevkenbohm et alj (119991) . 
iRreykenbo iim et all (2004), and Fritz et al. (20061). Fig. ^ 
shows the residuals for fitting an absorbed f dcut model to 
the averaged broad band spectrum for the first part of orbit 1 
(lA in Fig. [Hand Table [U, also including a Gaussian line to 
model the narrow Fe K emission feature. The reduced for 
this fit is Xred~24 f^i" degrees of freedom (dofs) and an ab- 
sorption line is clearly visible in the residuals at ^ 25-30 keV, 
which we interpret as a CRSF To model the CRSF we used a 
line shape with a Gaussian optical depth profile which is the 
basis for past modeling of RXTE data (local XSPEC model 
gauabs). The line is multiplicative with respect to the con- 
tinuum leading to the basic spectral model 

/(£) = F(E) X exp (- -^e'^- ) , (3) 

V /TTCTcyc ^ 

where F(E) is given by Eq. |2] and r^yc, £cyc, and (J eye are 
the optical depth, centroid energy, and width of the cyclotron 



line. While including this component improves Xied ^'^ 3. 5 for 
86 dofs a strong residual structure remains around 13keV 
(Fig. [2]:). Adding a single Gaussian emission line at this en- 
ergy results in a good overall fit with Xred ~ 

1.13 for 83 dofs 

(Fig. |2ji). The same model was also successfully applied to 
describe the other three time-averaged spectra with the re- 
sulting best fit parameters listed in Table [T] We note that 
the line parameters do not significantly depend on the con- 
tinuum model chosen: fits using a negative-positive exponen- 
tial continuum (local XSPEC model npex, Mihara 1995) in- 
stead of the power law with Fermi-Dirac cutoff lead to the 
same cyclotron line parameters. Compared to a NPEX or 
highecut continuum we achieved slightly better spectral 
fits using a power law with Fermi-Dirac cutoff. BOO used 
blackbody emission to describe the spectrum between 0.1 
and 2keV. Our spectral data does not extend below 3.5 keV 
and we d o not detect an influenc e of a black-body at lower 
energies. ' Becker & WoU dlOOTh are using a Comptonized 
Bremsstrahlung spectrum to calculate the continuum. They 
state that the contribution of blackbody emission is negligible 
in their spectrum. 

To better understand the observed 13keV feature, we 
checked for correlations for all possible pairs of spectral pa- 
rameters using X' contours. We found no correlation between 
the parameters characterizing the broad 13keV feature and 
the cyclotron line parameters or the roll-over of the continuum 
model. As expected, however, a dependence on the power law 
index F (and therefore also on A^h) was clearly detectable. An 
increasing F or A^h value resulted in a lower centroid energy 
of the 13 keV feature. We conclude that this feature is part of a 
more complex continuum rather than the simple phenomeno- 
logical approach used here. 

A broad feature at energies below the fundamental cy- 
clotron line, mostly around ^ lOkeV, has also been ob- 
served in several other cyclotron line sources, e.g., with 
RXTE in MXB 0656-072 (McBride et al. 20 0^), Her X- 
l,4U 1626-67, 4U 1907+09, and 4U 1538-52 (ICoburn et alJ 
I2OO2I) . and with Ginga in 4U1538-52, 4U 1907-09, and 
V 033 1+53 (lMiharall995h. A weak lOkeV ab sorption feature 
was also reported bv ISantangelo et aTl (119981) for a different 
BeppoSAX observation of Cen X-3. Such residuals are usu- 
ally best described by an additional broad absorption line to 
improve the fits. In some cases a broad emission line can de- 
scribe the data better. In our data we find that an emission line 
at 13 keV describes the observed feature much better. Substi- 
tuting the 13keV emission line by a Gaussian absorption line 
at ^ 8keV leads to Xred - ^--^l compared to 1.13 for 83 
dofs . 

As described in |T] BOO observed quasi-simultaneously 
with BeppoSAX (11 hours overlap with observation 1 A from 
Table[T]i. A phase averaged analysis was carried out by these 
authors using a power law with high-energy cutoff, smoothed 
around the cutoff energy for the continuum. The CRSF pa- 
rameters, with Ecyc = 30.6 ±0.6 keV and Ccyc = 5.9 ±0.7keV, 
are comparable to those presented here . Their measured 
value of A^H = 1.95 ± 0.03 x 10^^ cm"^ is still within the er- 
ror bars of our RXTE observation of A^h = 1 .6^q jcm"^. Their 
photon index of F = 1.208 ± 0.007 is much higher than our 
value of F = 0.92;!;q q4 due to our model definitions. We ap- 
plied the BOO model to our data and obtained consistent re- 
sults within error bars. Due to a known deviation in the 
slope of the F-ind ex between RXTE -PCA and BeppoSAX 
dKirsch et al.ll2005h . our resulting value of F = 1 .26 ± 0.03 for 
the BOO model is marginally consistent. 
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4.2. Physical model for CRSF 

Ideally, models based on a physical description of the 
accretion column above the magnetic poles of the neu- 
tron star should be used to fit the shape of the cyclotron 
line and the spectral c ontinuum. Despite recent successes 
dBecker & Wolffil2007l) . however, a fully self-consistent so- 
lution to the problem of continuum formation and radiative 
transfer in neutron star accretion columns is not ye t avail- 
able . Based on earlier work bv jArava & Hardin d ( Il999h 
and lArava-Gochez & Hardind (|2000|) . we have recently de- 
veloped a model for the self -consistent dete rmination of cy- 
clotron line profiles (Schonhe rr et al.ll2007h based on a pre- 
scribed continuum. Making use of a Monte Carlo approach 
and using the correct quant um-electrodynamical cross sec- 
tions, ISchonherr et al.l ( l2007h calculated the Green's function 
for the radiative transport for several possible geometries of 
the accretion column. By convolving an assumed input con- 
tinuum spectrum with these Green's functions, it is possible 
to calculate the photon spectrum emerging from the accretion 
column. The magnetic field is kept constant throughout the 
accretion column. This approach has been implemented into 
XSPEC, allowing the direct comparison of the model with 
observational data. We are presently working on an improved 
version of the model with magnetic field gradients. 

Here, we will present results from our analysis of the time 
averaged 20-70 keV HEXTE data with this model. We limit 
ourselves to the Green 's functions in the 1-0 geometry for 
the line forming region (Schonherr et al. l2007t iFreeman et alJ 
[1999; Isenberg, Lamb & Wang 1998), where the accretion 
column is described as a slab with the slab normal being par- 
allel to the magnetic field and where the continuum spectrum 
is illuminating the bottom of the slab. This setup mimicks a 
plane-parallel and thin emission region close to the surface. 
The continuum is again taken to be a power law modified by a 
Fermi-Dirac cutoff (Eq.|2]). We decided to excluded A^h from 
the continuum because it is neglectable in the chosen energy 
range A gravitational redshift z = 0.3 is assumed, a typical 
value for neutron stars (see A detailed analysis of several 
different sources, including phase resolved analysis, will be 
presented in a future paper 

Fig. [3^ shows the residuals after modeling the HEXTE 
spectrum with the continuum model without taking the cy- 
clotron line formation into account. Introducing the cy- 
clotron line model and refitting significantly reduces the 
residuals (Fig. [3J5), although some s tructure remains. As 
has been shown by ISchonherr et al.l (12007) . this structure 
is caused by the fact that the predicted self-consistent cy- 
clotron line shapes are generally too deep and in addition 
have significant emission wings, assuming the simplified 
scenario of a constant magnetic field and simple geometry. 
These wings are due to so-called spawned photons, which 
are the result of the stepwise de-excitation of electrons from 
higher L andau le vels after their excitation through a scattering 
photo n dArava & Harding 1999; Arava-Gochez & Harding 
12000). One way of reducing the line depth and strength of the 
wings is by invoking a partial covering approach, where only 
part of the seed photon spectrum is seen through the accretion 
column, while the remainder is observed without being mod- 
ified by the column. Note that this approach is not unique. 
Other geometries, such as a seed photon source distributed 
throughout the accretion column, a column with a tempera- 
ture gradient, or a vertically varying magnetic field, could also 
result in shallower cyclotron line profiles with less dominant 



10 r: 
5 r 
Or 
-5 r 
-10 r 



X 



4 
2 

?^ 



-4 



+ 



, + ++ +++++ 

++* ''''%V" 



b) 



rc) 



20 50 
Channel Energy [keV] 

Fig. 3. — Residuals for fitting data with CRSF shapes obtained for Monte 
Carlo simulations applied on phase averaged HEXTE data from 20-80 keV. 
(a) continuum f dcut*powerlaw with Xj^j = 24.0 (45 dofs). (b) included 
cyclomc model for Slab 1-0 geometry showing emission wings (Xred ~ 
with 41 dofs). (c) added partial covering to reduce emission wings (x^d = 1-4 
with 39 dofs). Note that the residuals in (a) are sUghtly different than those 
shown in Fig.[2]a since only the HEXTE data were modeled. 



wings (ISchonherr et alj|2007h . 

The result of fitting the data with a partial covering model 
is shown in Fig. [3t- This model describes the data well 
(xL= 1-35 for 40 de grees of freedom). The best fit param- 
eters are B = 3.46^0 03 x 10^^ G for the magnetic field, kTe = 
T.04-^°il keV for the electron temperature, tj = 2.4![] ;|]^ x 10"^ 
for the Thomson optical depth of the continuum, not to be 
confused with the optical depth of the gauabs model. Note 
that the optical depth in the lines can be larger by a factor of 
more than lO'* (Arava & Harding 1999). The cosine of the 
angle, 9, between our line of sight and the magnetic field is 
cos(6l) = 0.94!;j{J°, i.e., 9 = 20!^" degrees. 

To check the consistency of our result with that determined 
in ^4.11 we use the best-fit parameters from abo ve to calcu- 
late th e centroid energy of the cyclotron line ( Schonherr et alj 
I2OO7I Eq. 3) 



■^cyc ■ 



: nieC 



1+2B/Bcritsin^(0)-1 
sin^(6i) 



1 

T+z 



(4) 



giving E^yc = 30.9![J:^keV for data set lA. The width of 
the CRSF c an be estimated from Doppler broadening (e.g., 
iMeszaros & Nagel l985) 



81n(2)fc7; 



|cos(6')|£'c: 



(5) 



In our case we obtain a width of 8.0;!^ 9 keV. Both results lie 
within the uncertainties of Table [T] The CRSF depth in this 
model is described by the Thomson optical depth and cannot 
be directly translated into the gauabs depth due to strong 
dependencies of the scattering cross section on the emission 
angle 9. 
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Nh [10^^ cm-'] 

Fig. 4. — contours for the four phase averaged spectra presented in 
Table[T]and the pre-ecHpse dip. For each data set the 1- and l-a confidence 
levels are shown. Note that the individual Nw-V correlation is a non-physical 
effect described in the text. 

4.3. Spectral evolution over the orbit 

This data set allows us to perform the most complete and 
detailed analysis of the orbital dependence of Cen X-3's broad 
band spectrum to date. We extracted PCA and HEXTE spec- 
tra from the same data modes as for the time averaged analysis 
(ED but using shorter GTIs with their duration mainly con- 
strained by observation gaps due to SAA passages and earth 
occultations. This selection procedure leads to spectra with a 
typical exposure time of ^ 3.5 ks (A(/)oi-b ^ 0.03). Data sub- 
sets with considerably longer observation times were further 
split into subsets of ^ 3.5 ks duration. We ignored data taken 
during the eclipse and we also ignored the data set at the end 
of the pre-eclipse dip in orbit 1, since no HEXTE data were 
available. The resulting broad band spectra were then mod- 
eled with t he sa me model used to describe the time averaged 
spectra in ^4.11 The chosen energy range was 3.5-23 keV for 
PCA and 18-60keV for HEXTE. 

Fig.[Tf'-e shows selected time-resolved spectral parameters. 
As already indicated by the time averaged results (Table [TJ, 
the two orbits show a somewhat different behavior with re- 
spect to the two luminosity levels observed in each case. Dur- 
ing the first orbit the parameters are essentially stable within 
the uncertainties over most of the orbit. A^h increases near 
eclipse ingress and egress during both orbits (Fig. [TJ») as ex- 
pected from a n increasing amount of ste llar material in the 
Une of sight (IClark. Minato & Mil [19881) . During the pre- 
eclipse dip in orbit 1, A^h is enhanced by a factor of ~ 2. F 
also increases during the pre-eclipse dip (Fig.[T];). Error con- 
tours for the four phase averaged observations from ^4.1l and 
for the pre-eclipse dip in the first orbit are shown in Fig. |4] 
The observed A^jj-F correlation for the individual data sets 
is real and can not be sufficiently explained by the calculated 
error contours. This correlation is of unknown origin. The 
A^f|-F correlation observed for each individual data set is an 
artificial effect from the continuum modeling. When model- 
ing the spectrum one achieves similar values for different 
A^f|-F combinations, since increasing A^jj can compensate for 
a softer power law component. Therefore a non-physical cor- 
relation in A^jj and F is expected. 

During the second orbit we observe a change of the gen- 
eral behavior of many of the spectral parameters that is ap- 
parently related to the mid-orbit luminosity drop. A^h stays 
relative constant prior to mid-orbit and rises continuously af- 
terward (Fig. ^i>). Analysis with a frozen F value to take the 
A^H-r-correlation into account, showed the same tendency in 
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Fig. 5. — (a) values of Nw as a function of orbital phase (crosses) and cal- 
culated model for parameters mentioned in text (dashed line), (b) Difference 
of A^H and wind model. 

A^H during the second orbit. As in the pre-eclipse dip of or- 
bit 1, the power law F slightly increases with A^h (Fig. [It). 
The CRSF parameters variability is slightly more pronounced 
during the second orbit (Fig. [TJ/ and Fig. [D). Note that the 
cutoff and folding energies do not show any significant change 
throughout the whole observation. 

4.4. 'Wind model 

We calculated a simp le radiation driven wind model for A^h 
based on the model by ICastor. Abbott & KleinI (I1975L here- 
after CAK), although we note that this assumption of a radia- 
tion driven wind is only an approximation, as it is very likely 
that the wind is disrupted by the strong X-rays present in the 
system dWoidowski. Liedahl & Sakdl2001l: IWoidowski et alj 
I2003L and references therein). We assumed a mass loss, 
m, for the donor star and the CAK velocity profile y(r) = 
Voo(l -R*/r)^ , where Voo is the terminal velocity. This gives 



nnir) = 



Ankr^vir) 



(6) 



for the radial particle density profile, where A: is a conversion 
factor from mass to effective proton particle density, assuming 
a wind with cosmic abundances, a system inclination of 90°, 
and ignoring the orbital eccentricity (e < 0.016). We then 
integrate over the density along the line of sight for a specific 
position in the orbit: 



A^H 



nn{l')dl\ 



(7) 



where /ns = Rdfio^{4>oib) is the position of the neutron star, pro- 
jected to the celestial plane and where = 19 Rq = L58_^ 
is the separation between neutron star an d donor (iBlondin 
19941) . Assuming a typical value of (3 = 0.8 (iFriend & Abbott 
19861) we obtain good results for m = 7 x 10~^MQyr"' and 
Voo = 1200kms~ (Fig. [Sjj). For better comparison we also 
show the difference between data and model (Fig. [Sj?). Due 
to the definition of the model a strong correlation between in- 
creasing mass loss and terminal velocity is present, such that 
absolute values cannot be determined by this approach. 

Both orbits show a different behavior in the evolution of 
A^H- During the first orbit, where A^h is mostly constant, we 
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Fig. 6. — Normalized PCA and HEXTE pulse profiles for four energy 
bands: (a) 4.8-10.2keV, (b) 10.2-13.1 keV, (c) 13.1-18.2keV and (d) 18.3- 
49.3 keV. A secondary peak at pulse phase ~ 0.7 is visible for lower energies. 

see a symmetrical deviation during eclipse ingress and egress, 
when our simple approach fails to model the data. A similar 
deviation is seen at the beginning of the second orbit. After 
mid-orbit A^h keeps increasing for the remainder of the orbit, 
indicating an increased absorption in the seco nd half of the 
orbit. A similar behavior has been simulated bv lBlondin et alj 
([1990.) and probably indicates the presence of an accretion 
wake in this system. A small bump around (porh = 0.3-0.4 can 
be observed in both orbits, more pronounced in the second 
one. This bump also exists in the simulation and is another 
indication of a tidal wake. 

5. PHASE RESOLVED SPECTROSCOPY 
5.1. Pulse profile 

Pulse profiles with 32 phase bins for 4 energy bands were 
extracted for each of the 12 Obslds. We decided to use locally 
determined pulse periods for each ObsID and the ephemeris 
from Nagase et al. ( 1992), after using the ephemeris and pe- 
riod from BOO did not result in consistent pulse profiles. All 
pulse profiles showed a similar shape with a main peak at 
pulse phase (0puise) ^ 0.3, which allowed us to combine them 
after phase alignment. Figure |6] shows the resulting pulse 
profiles for four different energy ranges. A second weak 
peak a t pulse phase of ~ 0.7 is visible at energies below 
18 keV. 'N agase et alJ (119921) observed Cen X-3 in 1989 March 
with Ginga when the source had a 1 - 37 ke V luminosity of 
5 X lO"'^ ergs"\ a factor of two fainter than during this RXTE 
observation, and presents pulse profiles for 9 energy bands. 
At lower energies they observe a double peaked profile, where 
both peaks have equal intensities below 7 keV. Above 18 keV 
the pr ofile is comparable to our results with only one main 
peak. lAudlev et atl (12001 !) showed that the 2-25 ke V pulse 
profile is highly variable, changing between an asymmetric 
double peaked, as in our observation, and a much more com- 
plex pulse profile with multiple peaks, depending on the count 
rate of the binary system. 

We used the pulse profiles presented in Fig.|6]to create hard- 
ness ratios for each energy band using the lowest energy range 
as baseline (Fig. [TJ. At <j>puhe = 0.7 we see the decrease in 
hardness of the secondary peak toward higher energies. The 
rise and decay of the main peak (0puise = 0.1 and 0.5) show 




Fig. 7. — Normalized hardness ratios. Eachpanel has been calculated 
by taking the according pulse profile from Fig.|6]divided by panel (a) from 
Fig.|6] The pulse profile is indicated in the bottom panel (dashed line). Top 
right corner show the energy ranges. 

peaks in the hardness ratio, a clear indication that the lower 
energy emission region is beamed more sharply. It is not clear 
if a third peak, visible at (ppuhe = 1 -0, originates from a real 
physical feature or is part of an artificial wing from a phase 
continuity in the plot. 

5.2. Spectral evolution over the pulse profile 

Extraction of phase resolved spectra for each individual Ob- 
sld with 32 phase bins resulted in large uncertainties for the 
individual parameters. We therefore created phase resolved 
spectra for each ObsIDs with 8 phase bins and observed sim- 
ilarities in the behavior of the spectral parameters, e.g., an in- 
crease of the CRSF centroid energy during the rise of the main 
peak. We decided to sum all ObsIDs and use 32 phase bins 
resulting in an integration time of ~12ks for a phase analy- 
sis. In the pulse minimum we binned the phases by a factor 
of 2, and even by a factor of 6 for the secondary peak, to be 
able to constrain the CRSF feature. This results in a total of 
24 phase bins on which we performed a spectral analysis us- 
ing the model defined in ^4.11 Fig.|8]shows the best-fit values 
of selected parameters over the pulse period. Compared to the 
phase averaged analysis, not only Nh and F but also the CRSF 
parameters change drastically over the pulse period. 

The pulse profile can be roughly divided into three regions 
where the continuum parameters are significantly different: 
main peak, secondary peak, and pulse minimum. The col- 
umn density A^h is constant with a value of 1 .5"^^ x 10^^ cm"^ 
during the main peak. Between main and secondary peak, 
this value suddenly jumps up, to 2.2;|;q3 x 10^^ cm"^. con- 
tours indicate that this increase is not due to statistical fluctua- 
tions. During the pulse minimum the A^h value reaches a min 

r+0.5 . 



imumof0.5!!J:5 10 



cm"^. Another interpretation would be 
that A^H follows the shape of the pulse profile and is some- 
how capped during the main peak. We have not identified a 
physical explanation for this observational fact which might 
be caused by the limitations of the current continuum model- 
ing. The power law index F also changes strongly over the 
pulse phase. It increases during the rise of the pulse up to 
a value of I-O^q [ and then decreases to 0.5;!;q3 for the dura- 
tion of the main peak. After the main peak the value jumps 
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Fig. 8. — Selected spectral parameters over the pulse profile (dotted line). 
The last phase bins are binned. Error bars indicate Icr deviations. 



back up to l-2;!;o j. The CRSF centroid energy also shows 
a drastic change throughout the pulse. During the rise of 
the main peak the value increases from ^ 30keV to almost 
40keV, slowly decreasing throughout the main peak and de- 
scent below the phase average value of less than 30 ke V. The 
evolution of the CRSF throughout the pulse is similar to the 
shape of the pulse profile with a phase shift of A^puise = 0.2 
(ignoring the two outliers in the pulse rise we still obtain 
A^puise = 0.1). The width (Tcyc of the CRSF approximately 
follows the evolution of the centroid energy over the main 
peak, and drops to during the off-pulse, an indication 
that the CRSF is not well constrained. The depth Tcyc is 



slowly decreasing throughout the whole pulse and is also 
much less constrained during the pulse minimum. Compared 
to the other known sources featuring a variation of the CRSF, 
e.g., HerX- 1 (Gruber et al. 2001), Vela X-1 (La Barbera et aD 
^03"; Ki-evkenbohmet al. 2002), 4U 0352+309 (Coburn.eLai] 
2001), and GX 301-2 ( Kreykenbohm et al. 2004), Cen X-3 is 



the only known source to show an increase of the centroid 
energy during the ascent of the main pulse. 

BOO also provided a phase resolved analysis, but they di- 
vided the pulse into only 4 parts: ascent, maximum, descent of 
the main peak, and minimum, including the secondary peak. 
They also observe an increase in the CRSF centroid energy up 
to 36.6;I;2 4 keV in the ascent. We combined the phase resolved 
spectra according to the phase bins in BOO and applied their 
model on our results. We can confirm the results from BOO 
within the error bars with small deviations for F. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we presented the data analysis of the HMXB 
Cen X-3 observed for two consecutive orbits. We conducted 
phase averaged analysis for different parts of the orbit, distin- 
guished by a drop in the overall count rate. A decrease in the 
2-lOkeV luminosity between first and second orbit by 20% 
can be observed and could be caused by a decrease in the over- 
all amount of matter accreted onto the neutron star. In addi- 
tion, pulse phase resolved spectroscopy was performed. Here, 
the long duration of the observations allowed us to study the 
pulse phase dependence of the cyclotron line parameters in 
unprecedented detail. In the following two sections, we sum- 
marize the major results of these two studies, starting with the 
orbital variability of the X-ray spectrum. 

6.1. Orbital Variability of the X-ray Spectrum 

The analysis for individual satellite orbits shows the evo- 
lution of the spectral parameters throughout both orbits. We 
see a difference between both orbits, not only in the lumi- 
nosity, but also in the evolution of A^h and F, both correlated 
to each other Comparison with a simple CAK wind model 
showed that, while reproducing the overall trends in A^h, such 
a model is not sufficient to describe the observed data fully, 
regardless of inaccuracies during e clipse ingress and egress. 
IWoidowski. Liedahl & Sakol (1200 lb use a similar approach on 
a Cen X-3 data set from ASCA with significantly different re- 
sults. Probing the parameter space they find best fit values for 
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with S = 1.56. Applying these parameters to our model we 
find a minimum A^H value of 3.2 x lO^^cm"^, 1.5 times higher 
than our observed data. Also the change of A^h during ingress 
and egress of the eclipse is not as drastic as observed here. 
With a fixed /3 = 0.8 and S = 0.7 we obtain the results shown 
in Fig.|5] 

In our data set a small increase in A^h is seen in the first 
half ((f)oTb = 0. 3-0.4) in both orbits, more pronou nced in the 
second orbit. iBlondin. Stevens & Kalknanl (Il991h used two- 
dimensional numerical simulations based on the CAK model 
to calculate the column density in massive X-ray binary sys- 
tems with a tidal stream and observes such an increase in A^h 
around 0orb = 0.4-0.5. They interpret this as a leading bow 
shock in front of the accretion stream which is extremely sen- 
sitive to orbital parameters, e.g., distance between neutron star 



10 



Suchy et al. 



and donor or the intensity of the stellar wind. A possible ex- 
planation of our observed increase is that it corresponds to this 
bow shock observed in the sim ulations. I n the second half 
of the orbit, iBlondin. Stevens & Kallmi^ (Il991h see a gen- 
eral increase in A^h, comparable to our observations during 
the second orbit. This increase could be due to material in the 
line of sight, trailing the tidal wake. Taking into account that 
the continuum parameters do not change significantly over 
the orbit we can rule out that the system changed its state. 
Similar incre ases have been observed in othe r sources, e.g., 
4U 1700-37 dHaberl White & KallmanI [19891) and Vela X-1 
(iHaberl & W hite 1990h. both systems are believed to have a 
tidal streams. In these cases, the difference was more than an 
order of magnitude, whereas here the increase is on the order 
of a factor of 2-3. We point out that we also observe a change 
in A^H with a fixed photon index so that the A^e-r-correlation 
alone can not be responsible for this rise. The fact that this 
increase is only visible during the second orbit could imply 
that the overall luminosity has an influence on the visibility of 
the accretion wake. 

We do not rule out other possible expla nations for the dif- 
ferences in orbit 1 and 2. iBlondin et alj (Q990) showed in 
their simulations that the column density throughout the orbit 
changes between consecutive orbits due to variability in the 
accretion flow. The drop in luminosity of 20 % between 
orbit 1 and 2 indicates the decreas e of m. In the simulations 
IBlondin. Stevens & Kallmaiil (119911) observe instabilities for 
lower luminosities, i.e., less photon pressure from the donor 
star, trailing the tidal wake. These instabilities can create 
small pockets of denser material which also could be respon- 
sible for the observed increase in A^h- 

Another hint for the existence of multiple material clumps 
is seen in the light curve. Compared to the first orbit, the 
second orbit shows multiple pre-eclipse dips. In the first orbit, 
during the pre-eclipse dip A^h increases by a factor of ^ 2. 
Multiple dips close together could cause an overall increase 
of A^H throughout the second orbit. One small drawback to 
this theory is that the increase in A^h starts directly after the 
mid-orbit break, whereas the pre-eclipse dips are occurring at 
0oib > 0.7 of the second orbit. Therefore we do not think that 
pre-eclipse dips can be exclusively responsible for the rise. 

6.2. Spectral and Pulse Phase Variability 

We divided each orbit into two parts, before and after mid- 
orbit break, and extracted phase averaged spectra for each 
part. For the analysis we applied a power law with a Fermi- 
Dirac cutoff as continuum with an additional iron emission 
line. A CRSF at ^ 30keV is visible in all data sets, consis- 
tent with previous observations. The widely used local model 
gauabs gave us best fit results for this feature. No significant 
changes in the spectral parameters could be observed through- 
out the orbit. An additional broad Gaussian emission feature 
at 13 keV had to be introduced to get reasonable residuals. We 
checked for correlations of this feature with other model pa- 
rameters and determined that it is probably an artifact of the 
current continuum modeling. We conclude that the 13keV 
feature is equal to the 10 keV absorption feature that has been 
observed in multiple different sources and with different in- 
struments. Further investigation of this feature is necessary, 
but this is beyond the scope of this paper 

We also applied a newly de veloped physical model, 
cyclomc (ISchonherr et al.l l2007h . based on Monte-Carlo 
simulations to part of the phase averaged data (data set lA). 
This model is self-consistent in determining the CRSF param- 
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Fig. 9. — Calculated magnetic field strengths (crosses) and BOO results 
(triangles) using the CRSF centroid energy. The dotted line show the single 
dipole model using par ameters from BOO. The solid line show the best fit of 
the single dipol model to our data. 

eters using the Green's function for radiative transport. For 
all geometries of the line forming region available in the cur- 
rent implementation of the model, we generally find that the 
predicted CRSFs are too deep and possess significant emis- 
sion wings, leading to difficulties when matching the observed 
spectra. The slab 1-0 geometry, where the continuum is illu- 
minating the bottom of the accretion column, has the small- 
est emission wings of the geometries tested and is therefore 
a good candidate for a first approach. Assuming partial cov- 
ering of the continuum leads to a further reduction of the re- 
maining emission wings. This approach is motivated by the 
idea that only a part of the continuum contribution emerges 
from the line forming region. The results are encouraging and 
in good agreement with the previous phenomenological ap- 
proach. It is planned to improve this model by, e.g., adding 
magnetic field that are gradually changing over the accre- 
tion column, and to apply it to different sources with known 
CRSFs as part of a future paper. 

For the phase resolved analysis we created pulse profiles 
in different energy bands over the whole observation. They 
show a main peak and a much smaller secondary peak which 
is only seen at lower energies. Compared to other sources, 
e.g.. Vela X-1 or GX 301-2, the pulse profile is rather smooth 
and does not show a complex shape at lower energies. The 
main peak is more pronounced at lower energies and slightly 
asymmetric, steeper in the rise, at higher energies. Comparing 
with Nagase et al. (1992), where a clear double peaked profile 
at lower energies has been observed, we assume that this is a 
luminosity de pendent effect on the pulse profile. During the 
observation of iNagase et al.l (Il992h . Cen X-3 was fainter than 
during our second orbit and no changes in the pulse profiles 
have been observed throughout our observation. 

For each bin in the pulse profile we extracted a spectrum to 
study the evolution of the parameters over the pulse. For bet- 
ter statistics we combined some of the data in the pulse min- 
imum, resulting in 24 phase bins overall. The continuum pa- 
rameters show some significant changes throughout the pulse 
and the general A^e-r-correlation seems not to be valid during 
the main peak. 

The CRSF shows significant changes throughout the pulse. 



Cen X-3 over two binary orbits 
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The shape of the changes is very similar to the pulse profile 
shifted by A^puise = 0.1, not including the two maximum val- 
ues during the rise. These outliers have rather large uncer- 
tainties and should be handled with caution in the interpreta- 
tion. From the main peak on, E^yc slowly decreases to below 
30keV where it stays for the remainder of the pulse. Fig. |9] 
shows how the magnetic field deduced from these measure- 
ments varies as a function of pulse phase. The figure also 
shows the earlier BeppoSAX measurements of BOO, which 
are consistent with our much finer resolved data. 

BOO argued that the large variation in Z?-field with pulse 
phase could not be due to a variation of the height in the accre- 
tion column at which the observed radiation originates, since 
this variation would imply height differences of kilometers. 
BOO instead argue that, based on earlier studies of the pulse 
shape of Cen X-3 by Leahy ( 1991), the polar cap of Cen X-3 
is rather large. If this assumption is true, then it is likely that 
different locations on the neutron star's surface are observed 
over one rotation of the neutron star Using the assumption 
that the observed X-rays are only produced by one of the two 
magnetic poles, and furthermore assuming that the observed 
X-rays are all coming from a location at a constant geograph- 
ical colatitude of the neutron star (taken by BOO to be the 
system's inclination), BOO then show that the variation of the 
Z?-field at the neutron star's surface produced by a magnetic 
dipole offset by 0.1 /?ns from the rotational axis of the neu- 
tron star can explain the Z?-field variation seen by BeppoSAX. 
We reproduce their best-fit model in Fig. |9l While the model 
is sufficient to explain the BeppoSAX -data, the higher reso- 
lution i?XTE-data clearly show that it is not a viable explana- 
tion for the observed Z?-field variation. We extended the ap- 
proach of BOO by allowing for observations of both magnetic 
poles and by also removing the constraint that the colatitude 
at which the X-rays originate equals the inclination. Neither 
of these approaches resulted in statistically satisfying descrip- 
tions of the variation of B with </>puise ■ Even when relaxing the 
model assumptions further, by using a skew-symmetric mag- 
netic dipole, which approximates higher multipole moments 
for the B-field on the surface, and assuming that the X-rays 
come from two different locations on the neutron star, no sat- 
isfying description of the observed line variation was found. 
We therefore conclude that the model of BOO is not a viable 



description of the high resolution B-field variation seen here. 

We note that an alternative s tudy of the pulse profile of 
Cen X-3 has been presented bv iKraus et alj ([T996 b. Similar 
to the assumptions outlined above, these authors show that 
the shape and energy dependence of the pulse profile can be 
explained by emission from a distorted magnetic dipole. Tak- 
ing into accoun t the relativi s tic ph oton propagation close to 
the neutron star, iKraus et aTl d 19961) show that, in order to ex- 
plain the pulse profile, both magnetic poles must contribute 
to the observed emission. Modeling observations of Cen X-3 
made at a si milar lumino sity state as the one during the RXTE 
observation. iKraus et alJ 11996) show that these assumptions 
yield an offset of approximately 10° for both poles, and that 
both poles contribute equally to the observed flux. In these 
models, they assume that the polar caps radiation characteris- 
tic is a pen- plus pencil-beam pattern. A major result of these 
models is that in this decomposition the main pulse emission 
is dominated by the magnetic pole facing toward the observer, 
while the wings of the main pulse are dominated by the pole 
facing away from the observer Although they do not take 
the height of the accretion column into account, we believe it 
is likely that for geometric reasons the X-rays observed from 
the magnetic pole on the neutron star hemisphere facing away 
from the observer originate at a distance farther away from 
the neutron star surface, and thus in a region of lower Z?-field. 
Together with the possibility of slightly different surface mag- 
netic strengths at both poles, we therefore deem it likely that 
the observed B-field variation over the X-ray pulse is due to us 
observing the two accretion columns at different heights, al- 
though further modeling of such a geometric setup is clearly 
required before a final answer on the origin of the B-field vari- 
ation can be given. 
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